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The Giant Planets 



Tristan Guillot 
Observatoire de la Cote d'Azur 

Daniel Gautier 
Observatoire de Paris 

Abstract 

We review the interior structure and evolution of Jupiter, Saturn, Uranus and Neptune, 
and extrasolar giant planets with particular emphasis on constraining their global composition. 

KEYWORDS: Giant planets, extrasolar planets, Jupiter, Saturn, Uranus, Neptune, planet 
formation 



1 Introduction 

In our Solar System, four planets stand out for their sheer mass and size. Jupiter, Saturn, Uranus 
and Neptune indeed qualify as "giant planets" because they are larger than any terrestrial planet 
and much more massive than all other objects in the Solar System, except the Sun, put together 
(fig. [1]). Because of their gravitational might, they have played a key role in the formation of the 
Solar System, tossing around many objects in the system, preventing the formation of a planet in 
what is now the asteroid belt, and directly leading to the formation of the Kuiper belt and Oort 
cloud. They also retain some of the gas (in particular hydrogen and helium) that was present when 
the Sun and its planets formed and are thus key witnesses in the search for our origins. 

Because of a massive envelope mostly made of hydrogen helium, these planets are fluid, with 
no solid or liquid surface. In terms of structure and composition, they lie in between stars (gaseous 
and mostly made of hydrogen and helium) and smaller terrestrial planets (solid and liquid and 
mostly made of heavy elements) , with Jupiter and Saturn being closer to the former and Uranus 
and Neptune to the latter. 

The discovery of many extrasolar planets of masses from a few hundreds down to a few Earth 
masses and the possibility to characterize them by the measurement of their mass and size prompts 
a more general definition of giant planets. For this review, we will adopt the following: "a giant 
planet is a planet mostly made of hydrogen and helium and too light to ignite deuterium fusion" . 
This is purposedly relatively vague - depending on whether the inventory is performed by mass 
or by atom or molecule, Uranus and Neptune may be included or left out of the category -. Note 
that Uranus and Neptune are indeed relatively different in structure than Jupiter and Saturn and 
are generally referred to as "ice giants", due to an interior structure that is consistent with the 
presence of mostly "ices" (a mixture formed from the condensation in the protoplanetary disk of 
low refractivity materials such as H2O, CII4 and NH3, and brought to the high pressure conditions 
of planetary interiors -see below). 
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Terrestrial Jupiter Saturn Uranus Neptune Kuiper belt Oort cloud 



Figure 1: An inventory of hydrogen and helium and all other elements ("heavy elements") in the 
Solar System excluding the Sun (the Sun has a total mass of 332, 960 M©, including about 5000 
in heavy elements, 1 being the mass of the Earth) . The precise amount of heavy elements in 
Jupiter (10 - 40 M®) and Saturn (20 - 30 M®) is uncertain (see 

Globally, this definition encompasses a class of objects that have similar properties (in par- 
ticular, a low viscosity and a non- negligible compressibility) and inherited part of their material 
directly from the same reservoir as their parent star. These objects can thus be largely studied 
with the same tools, and their formation is linked to that of their parent star and the fate of the 
circumstellar gaseous disk present around the young star. 

We will hereafter present some of the key data concerning giant planets in the Solar System and 
outside. We will then present the theoretical basis for the study of their structure and evolution. 
On this basis, the constraints on their composition will be discussed and analyzed in terms of 
consequences for the models of planet formation. 

2 Observations and global properties 
2.1 Visual appearances 

In spite of its smallness, the sample of four giant planets in our Solar System exhibits a large 
variety of appearances, shapes, colors, variability... etc. As shown by fig. [21 all four giant planets 
are flattened by rotation and exhibit a more or less clear zonal wind pattern, but the color of 
their visible atmosphere is very different (this is due mostly to minor species in the high planetary 
atmosphere), their clouds have different compositions (ammonia for Jupiter and Saturn, methane 
for Uranus and Neptune) and depths, and their global meteorology (number of vortexes, long- 
lived anticyclones such as Jupiter's Great Red Spot, presence of planetary-scale storms, convective 
activity) is different from one planet to the next. 

We can presently only wonder about what is in store for us with extrasolar giant planets since 
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Figure 2: Photographs of Jupiter, Saturn, Uranus and Neptune. 



we cannot image them. But with orbital distances that can be as close as 0.02 AU, a variety of 
masses, sizes, and parent stars, we should expect to be surprised! 



The mass of our giant planets can be obtained with great accuracy from the observation of the 
motions of their natural satellites: 317.834, 95.161, 14.538 and 17.148 times the mass of the 
Earth (IM^ = 5.97369 x 10^^ g) for Jupiter, Saturn, Uranus and Neptune, respectively. More 
precise measurements of their gravity field can be obtained through the analysis of the trajectories 
of spacecrafts during flyby, especially when they come close to the planet and preferably in a 
near-polar orbit. The gravitational field thus measured departs from a purely spherical function 
due to the planets' rapid rotation. The measurements are generally expressed by expanding the 
components of the gravity field on Legendre polynomials Pi of progressively higher orders: 



where VcxtifjO) is the gravity field evaluated outside the planet at a distance r and colatitude 9, 
i?eq is the equatorial radius, and Ji are the gravitational moments. Because the giant planets are 
very close to hydrostatic equilibrium the coefficients of even order are the only ones that are not 
negligible. We will see how these gravitational moments, as listed in table [1] help us constrain the 
planets' interior density profiles. 

Table [T] also indicates the radii obtained with the greatest accuracy by radio-occultation ex- 
periments. An important consequence obtained is the fact that these planets have low densities, 
from 0.688 gcni^^ for Saturn to 1.64 gcm^^ for Neptune, to be compared with densities of 3.9 to 
5.5gcm~^ for the terrestrial planets. Considering the compression that strongly increases with 
mass, one is led naturally to the conclusion that these planets contain an important proportion 
of light materials including hydrogen and helium. It also implies that Uranus and Neptune which 
are less massive must contain a relatively larger proportion of heavy elements than Jupiter and 
Saturn. This may lead to a sub-classification between the hydrogen-helium giant planets Jupiter 
and Saturn, and the "ice giants" or "sub giants" Uranus and Neptune. 

The planets are also relatively fast rotators, with periods of ^ 10 hours for Jupiter and Saturn, 
and ~ 17 hours for Uranus and Neptune. The fact that this fast rotation visibly affects the figure 
(shape) of these planets is seen by the significant difference between the polar and equatorial radii. 
It also leads to gravitational moments that differ significantly from a nul value. However, it is 



2.2 Gravity fields 




(1) 
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Table 1: Characteristics of the gravity fields and radii 





Jupiter 


Saturn 


Uranus 


Neptune 


M X 10"^^ [kg] 


18.986112(15)=' 


5.684640(30)'" 


0.8683205(34)"= 


1.0243542(31)"* 


i?eq X 10"^ [m] 


7.1492(4)"= 


6.0268(4)' 


2.5559(4)s 


2.4766(15)'^ 


iipoi X 10"^ [m] 


6.6854(10)"= 
6.9894(6)'" 


5.4364(10)' 
5.8210(6)'" 


2.4973(20)*^ 


2.4342(30)8 


R X 10"'' [m] 


2.5364(10)' 


2.4625(20)' 


p X 10"^ [kgm"^] 


1.3275(4) 


0.6880(2) 


1.2704(15) 


1.6377(40) 
0.3539(10)"* 
-0.28(22)"* 


J2 X 10^ 


1.4697(1^ 


1.632425(27)'" 
-9.397(28)'" 
0.867(97)'" 


0.35160(32)"= 


Ja X lO-* 


-5.84(5)'' 


-0.354(41)"= 


Je X 10^ 


0.31(20)'' 






P„ X 10-* [s] 


3.57297(41)j 


3.83624(47) V'^ 


6.206(4)' 


5.800(20)"" 


1 


0.08923(5) 


0.15491(10) 


0.02951(5) 


0.02609(23) 


C/MRl^ 


0.258 


0.220 


0.230 


0.241 



The numbers in parentheses are the uncertainty in the last digits of the given value. The 

value of the gravitational co nstant used to calculate the masses of Jupiter and Saturn is G = 

6.67259 X 10~*^Nm^kg"^ (jCohen and Tavloil Il987l l. The values of the radu, density and 

gravitational moments correspond to the one bar pressure level (lbar= 10^ Pa). 

" CampbeU and Svnnott (1985) 

'" Jacobson et al. (2006=) 

"= Anderson et al. (1987) 

"* Tvler et al. (1989) 

"= Lindal et al. ('1981') 

^ Lindal et al. ('1985') 

s Lindal (1992a) 

*" From 4th order figure theory 

' (2iicq + -Rpoi)/3 (Clairaut's approximation) 

j Davies et al. (1986) 

This measurement from the Voyager era is now in question and values up to 38826 s have been 
proposed (see ij |2.3|l 
' Warwick ot al. ( 1986) 
°" ~Warwick et al. (1989) 



important to stress that there is no unique rotation frame for these fluid planets: atmospheric 
zonal winds imply that different latitude rotate at different velocities (see ii l2.6p . and the magnetic 
field provides an other rotation period. Because the latter is tied to the deeper levels of the planet, it 
is believed to be more relevant when interpreting the gravitational moments. The rotation periods 
listed in Table [1] hence correspond to that of the magnetic field. The case of Saturn appears to be 
complex and is discussed in the next section. 

2.3 Magnetic fields 

As the Earth, the Sun and Mercury, our four giant planets possess their own magnetic fields, as 
shown by the Voyager 2 measurements. The structures of these magnetic fields are very different 
from one planet to another and the dynamo mechanism that generates them is believed to be 



related to convection in their interior but is otherwise essentially unknown (see IStevensonl 1X983 . 
for a review). 

The magnetic field B is generally expressed in form of a development in spherical harmonics of 
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the scalar potential W, such that B = 



-VW: 



W 



(7) 



n+1 



E 

rn=0 



{5^ cos(™0) + /i^T sin(m0)} P™(cos0). 



(2) 



r is the distance to the planet's center, a its radius, 9 the colatitude, (p the longitude and P,™ 
the associated Legendre polynomials. The coefficients g,™ and /i™ are the magnetic moments that 
characterize the field. They are expressed in magnetic field units. 

One can show that the first coefficients of relation ^ (for n — and n = 1) correspond to the 
potential of a magnetic dipole such that W — 'M. ■ r/r^ of moment: 



M ^a-" 



9lY 



1/2 



(3) 



Jupiter and Saturn have magnetic fields of essentially dipolar nature, of axis close to the rotation 
axis {gl is much larger than the other harmonics); Uranus and Neptune have magnetic fields that 
are intrinsically much more complex. To provide an idea of the intensity of the magnetic fields, the 
value of the dipolar moments for the four planets are 4.27 Gauss Rj , 0.21 Gauss Rp, . 23 Gauss Rjt , 



0.133 Gauss R^, respectively (jConnernev et all . Il982t lAcuna et all Il983t iNess et al.l Il986l Il989[ ) 



A true surprise from Voyager that has been confirmed by the Cassini-Huygens mission is that 
Saturn's magnetic field is axisymetric to the limit of the measurement accuracy: Saturn's magnetic 
and rotation axes are perfectly aligned. Voyager measurements indicated nevertheless a clear 
signature in the radio signal at lOh 39min 22s believed to be a consequence of the rotation of 
the magnetic field. Determinations of a magnetic anomaly and new measurements by Cassini 
have since considerably blurred the picture, and the interpretation of the measurements have 
become unclear , with a deep r otation period evaluated betw een that of Voyager and as slow as 
lOh 47min 6s (IGurnett et al. . 2005[ Giampieri et al. . 20061 ) (see also Galopeau and Lecacheux . 



2OOOI : Cecconi and Zarka . 2005 ). Note that models discussed hereafter have not yet included this 



additional uncertainty. 



2.4 Atmospheric compositions 

In fluid planets, the distinction between the atmosphere and the interior is not obvious. We 
name "atmosphere" the part of the planet which can directly exchange radiation with the exterior 
environment. This is also the part which is accessible by remote sensing. It is important to note 
that the continuity between the atmosphere and the interior does not guarantee that compositions 
measured in the atmosphere can be extrapolated to the deep interior, even in a fully convective 
enviro i iment: Processes such as ph ase separations (e.g. ISalpetei . 1973t Stevenson and Salpetei . 
I1977U iFortnev and Hubbardl. l2003f) . phase transitions (e.g. iHubbardl . Il989f ). chemical reactions 
(e.g. Feglev and Lodders . 1994( l can occur and decouple the surface and interior co mpositions . 
Furth ermore, imperfect mixing may also occur, depending on the initial conditions (e.g. StevensonL 
Il985h . 

The conventional wisdom is however that these processes are limited to certain species (e.g. 
helium) or that they have a relatively small impact on the global abundances, so that the hydrogen- 
helium envelopes may be considered relatively uniform, from the perspective of the global abun- 
dance in heavy elements. We first discuss measurements made in the atmosphere before inferring 
interior compositions from interior and evolution models. 
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2.4.1 Hydrogen and helium 



The most important components of the atmospheres of our giant planets are also among the most 
difficult to detect: H2 and He have a zero dipolar moment and hence absorb very inefficiently 
visible and infrared light. Their infrared absorption beco mes important only a t high pressures 



when collision- induced absorption becomes significant (e.g. iBorvsow et al.l . 119971 ). On the other 
hand, lines due to electronic transitions correspond to very high altitudes in the atmosphere, and 
bear little information on the structure of the deeper levels. The only robust result concerning 
the abundance of heliu m in a giant planet is by in situ measurement by the Galileo probe in the 
atmosphere of Jupiter ( von Zahn et al. . 1998h . The helium mole fraction [i.e. number of helium 
atoms over the total number of species in a given volume) is qhc = 0.1359 ± 0.0027. The helium 
mass mixing ratio Y {i.e. mass of helium atoms over total mass) is constrained by its ratio over 
hydrogen, X: Y/{X -\-Y) = 0.238 ± 0.05. This ratio is by coincidence that found in the Sun's 
atmosphere, but because of helium sedimentation in the Sun's radiative zone, it was larger i n 



the protosolar nebula: K 



pioto 



0.275 ± 0.01 and {X + Y) 



proto 



0.98 (e.g. iBahcall et all . 119951 ) 



Less helium is therefore found in the atmosphere of Jupiter than inferred to be present when the 
planet formed. We will discuss the consequences of this measurement later: let us mention that 
the explanation invokes helium settling due to a phase separation in the interiors of massive and 
cold giant planets. 

Helium is also found to be depleted compared to the protosolar value in Saturn's atmo- 
sphere. However, in this case the analysis is complicated by the fact that Voyager radio oc- 
cultations apparently l ed to a wrong value. The current adopted value is now Y = 0.18 — 0.25 
Conrath and Gautier ^ 2000f). in agreem ent with values predicted by interior and evolution models 



Guillotl . Il999bl : iHubbard et al.l . 1 19991) . Finally, Uranus and Neptune are found to have near- 



protosolar helium mixing ratios, but with considerable uncertainty. 



2.4.2 Heavy elements 



The abundance of other elements than hydrogen and helium (that we will call hereafter "heavy el- 
ements" ) bears crucial information for the understanding of the processes that led to the formation 
of these planets. 

The most abundant heavy elements in the envelopes of our four giant planets are O, C, N, S. It 
is possible to model the chemistry of ga ses in the tropospher e s from the top of the convective zone 
down to the 2000 K temperature level ( Feglev and Lodder^ . 1994 ). Models conclude that, what- 
ever the initial composition in these elements of planetesimals which collapsed with hydrogen onto 
Jupiter and Saturn cores during the last phase of the planetary formation, C in the upper tropo- 
spheres of giant planets is mainly in the form of gaseous CH4, N in the form of NH3, S in the form 
of H2S, and O in the form of H2O. All these gases, but methane in Jupiter and Saturn, condense 
in the upper troposphere, but vaporize at deeper levels when the temperature increases. Interest- 
ingly enough, noble gases are not expected to condense even at the cold tropopause temperatures 
of Uranus and Neptune. 

The mass spectrometer aboard the Galileo atmospheric probe has performed in situ measure- 
ments of Ar, Kr, Xe, GH4, NH3, H2S, an d H2O in th e trop osphere of Jupiter. C, N, and S were 
found to be oversolar by a factor 3 to 4 ( Wong et al. . 200J), which was not unexpected because 
condensation of nebula gases results in enriching icy grains and planetesimals. The surprise came 
from Ar, Kr, Xe, which were expected to be solar because t hey are diff i cult t o condense, but 
turned out to be oversolar by factors 2 to 4 (jOwen et al. . I1999I: IWong et al. I. I2004I) . One exception 
among these enriched species was neon, which was found to be significantly undersolar, but was 
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Table 2: Main gaseous components of heavy elements measured in the tropo- 
sphere of giant planets 





Species 


Mixing ratio/H2 


References 


Comments 


Jupiter 


CH4 


(2.37 ± 0.37) X lO"-" 


Wone ot al. f2004) 


GPMS on Galileo(l) 




NH3 


(6.64 ± 2.34) X 10"^ 


Wong et al. (2004) 


idem 




H2S 


(8.9 ± 2.1) X 10"^ 


Wong et al. (2004) 


idem 




H2O 


(4.9 ± 1.6) X 10"" 


Wong et al. (2004) 


idem; region not well mixed 






(6.1 ± 1.2) X 10"" 


Atreva et al. (1999) 


idem 






(1.84 ± 0.37) X 10"^ 


Atreva et al. (1999) 


idem 






(4.9 ± 1.0) X 10"" 


Atreva et al. (1999) 


idem 


Saturn 


CH4 


(4.3 ± 1) X 10"^ 


Flasar et al. (2005) 


cms on Cassini (3) 




NH3 


(1 ± 1) X 10"* 


Briggs and Sackett (1989) 


Ground-based microwave (4) 




H2S 


(2.2 ± 0.3) X 10"* 


Briggs and Sackett (1989) 


idem 


Uranus 


CH4 


(3.3 ± 1.1) X 10"^ 


Gautier et al. (1995) 


Compilation from ground-based 
observations 

Ground-based microwave (4) 




H2S 


(1 ± 1) X 10"* 


Briggs and Sackett (1989) 


Neptune 


CH4 


(3.3 ± 1.1) X 10"^ 


Gautier et al. (1995) 


Compilation from ground-based 
observations 

Ground based microwave (5) 




H2S 


(7.5 ± 3.25) X 10"* 


de Pater et al. (1991) 




H2O 


7.7 X lO(-l) 


Lodders and Feglev (1994) 


Inferred from CO (6) 



(1) Galileo Probe Mass Spectrometer aboard the atmospheric probe in Jupiter 

(2) The signal stopped at the 22 bar levels prior to have reached a constant value. It is currently 
believed that the region where the probe made measurementds was atypically dry and that the 
bulk abundance of H2O in Jupiter has not been measured. 

(3) Composite Infra Red Spectrometer aboard the Cassini spacecraft 

(4) Ground based measurements of the microwave continuum. The result is somewhat uncertain 
due to the difficulty to precisely estimate opacities of absorbing species 

(5) Ground based microwave measurements. An oversolar II2S abundance is required to interpret 
the depletion of NII3 in the upper troposphere 

(6) Inferred from the microwave detection of CO in the troposphere. Note that the validity of the 



approach is questionned bv iBezard et al.l (|2002l ). A large amount of water seems to be present 
anyway in the deep atmosphere of Neptune. The case of Uranus is still uncertain because it is not 
known so far if CO is present in the troposphere of the planet. 



predi cted to be so because of a capture by the falling hehum droplets (jRoulston and Stevenson . 



19951). Another exception was water, but this molecule is affected by meteorological processes, and 



the probe was shown to have fallen into a dry region of Jupiter's atmosphere. 

Specifically, CH4/H2 has been found oversolar in the four giant planets: the C/H ratio corre- 
sponding to the measured abundances is always higher than the solar C/H ratio, and in fact appears 
to be increasing with distance to the Sun. C/H is 3, 7.5, 45 and 45 times solar, in Jupiter, Saturn, 
Uranus and Neptune, respectively. Note that the quoted enrichments are subject to changes when 
the solar abundances tables are revised, which happens surprisingly frequently. 

Except for Jupiter, the determination of the NH3 abundance is more uncertain than that of CH4 
because it is model dependant. It is derived from fitting microwave spectra of giant planets which 
exhibit a continuum opacity, more difficult to model than absorption spectral lines. However, the 
N/H enrichment seems to be, so far, fairly constant from a planet to another, around a factor 2, 
so that C/N is higher in Saturn than in Jupiter, and still higher in Uranus and Neptune. 

H2S has been measured in situ in Jupiter, but in the three other giant planets its large abun- 
dance is derived from the requirement to deplete NH3 at deeper levels than the saturation one; 
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This scenario has been proposed a long time ago by iGulkis et al. ( 19781 ). It iniphes that S/H is 
substantiaUy oversolar in Uranus and in Neptune. 

H2O is difficult to measure in all four giant planets because of its condensation relatively deep. 
It was hoped that the Galileo probe would provide a measurement of its deep abundance, but the 
probe fell into one of Jupiter's 5-micro ns hot spot, what is now believ ed to be a dry region mostly 
governed by downwelling motions (e.g. Showman and Ingersolll . 1998[ ). As a result, and although 
the probe provided measurements down to 22 bars, well below water's canonical 5 bar cloud base, 
it is believed that this measurement of a water abundance equal to a fraction of the solar value is 
only a lower limit. 



2.5 Energy balance and atmospheric temperature profiles 

Jupiter, Saturn and Neptune are observed to emit significantly more energy than they receive from 
the Sun (see Table [3|). The case of Uranus is less clear. Its intrinsic heat flux Fi^t is significantly 
smaller than that of the othe r giant planets. Deta i led m odeling of its atmosphere however indicate 



that i^int ^ 60 erg cm ^s ^ ( Marlev and McKay . 19991 ). With th is caveat, all fou r giant planets 



can be said to emit more energy than they receive from the Sun. iHubbard (|l968l ) showed in the 
case of Jupiter that this can be explained simply by the progressive contraction and cooling of the 
planets. 



Table 3: Energy balance as determined from Voyager IRIS data'' 



Jupiter 



Saturn 



Uranus 



Neptune 



-21 



Absorbed power [10" J s"^] 50.14±2 
Emitted power [10^ Js"^] 
Intrinsic power [10^'' 
Intrinsic flux [J s^^ m 
Bond albedo [] 
Efi'ective temperature [K] 
1-bar temperature'' [K] 



11.14±0.50 0.526±0.037 0.204±0.019 



83.65±0.84 
33.5±2.6 
5.44±0.43 
0.343±0.032 
124.4±0.3 
165±5 



19.77±0.32 
8.63±0.60 
2.01±0.14 

0.342±0.030 
95.0±0.4 
135±5 



0.560±0.011 

034+" 

U.UOI _o.034 

042 

U.UI^ -0.042 

0.300±0.049 
59.1±0.3 

76±2 



0.534±0.029 
0.330±0.035 
0.433±0.046 
0.290±0.067 
59.3±0.8 
72±2 



After Pearl and Conrath 
^ Lindal (1992b) 



A crucial consequence of the presence of an intrinsic heat flux is that it requires high internal 
temperatures 10, 000 K or more), and that c onsequently the giant planets are fluid (not solid) 
(|Hubbard (1968,); see also lHubbard et al.l (|l995l )). Another consequence is that they are essentially 
convective, and that their interior temperature profile are close to adiabats. We will come back to 
this in more detail. 

The deep atmospheres (more accurately tropospheres) of the four giant p lanets are inde ed ob- 
served to be close to adiabats, a result first obtained by spectroscopic models (jTraftonl . 119671 ) . then 
verified by radio-occultation experiments by the Voyager spacecrafts, and by the in situ measure- 
ment from the Galileo probe (fig. [3]). The temperature profiles show a temperature minimum, in a 
region near 0.2 bar called the tropopause. At higher altitudes, in the stratosphere, the temperature 
gradient is negative (increasing with decreasing pressure). In the regions that we will be mostly 
concerned with, in the troposphere and in the deeper interior, the temperature always increases 
with depth. It can be noticed that the slope of the temperature profile in fig [3] becomes almost 
constant when the atmosphere becomes convective, at pressures of a few tens of bars, in the four 
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giant planets. 
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Figure 3: Atmospheric temperatures as a function of pressure for J upiter, Saturn , Uranus and 
Neptune, as obtained from Voyager radio-occultation experiments (see lLindall . ll992bl ). The dotted 



hne corresponds to th e temperature pr ofile retrieved by the Galileo probe, down to 22 bar and a 



temperature of 428 K (jSeiff et al.l . Il998f ). 



It should be noted that the 1 bar temperatures listed in table [3] and the profiles shown in Fig. [3] 
are retrieved from radio-occultation measurements using a helium to hydrogen ratio which, at least 
in the case of Jupiter and Saturn, was shown to be incorrect. The new values o f Y are f ound t o 
lead to increased temperatures by 5K in Jupiter and 10 K in Saturri (see i Guillot . 1999al ). 
However, the Galileo probe found a 1 bar temperature of 166 K ( Seiff et al1 . ll998h . and generally 
a good agreement with the Voyager radio-occultation profile with the wrong He/H2 value. 

When studied at low spatial resolution, it is found that all four giant planets, in spite of their 
inhomogeneous appearances, have a rather unifor m brightness temperat ure, with pole-to-equator 



latitudinal variations limited to a few kelvins {e.g. IngersoU et al. . 1995f) . However, in the case of 
Jupiter, some small regions are known to be very different from the average of the planet. This 
is the case of hot spots, which cover about 1% of the surface of the planet at any given time, but 
contribute to most of the emitted flux at 5 microns, due to their dryness (absence of water vapor) 
and their temperature brightness which can, at this wavelength, peak to 260 K. 



2.6 Atmospheric dynamics: winds and weather 

The atmospheres of all giant planets are evidently complex and turbulent in nature. This can for 
example be seen from the mean zonal win ds (inferred from clou d tracking), which are very rapidly 
varying functions of the latitude (see e.g. Ilngersoll et al. . 1995f ): while some of the regions rotate 
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at the same speed as the interior magnetic field (in the so-called "system III" reference frame), 
most of the atmospheres do not. Jupiter and Saturn both have superrotating equators (+100 
and +400 ms""'^ in system III, for Jupiter and Saturn, respectively), Uranus and Neptune have 
subrotating equators, and superrotating high latitude jets. Neptune, which receives the smallest 
amount of energy from the Sun has the largest peak-to-peak latitudinal variations in wind velocity: 
about 600 ms~^. It can be noted that, contrary to the case of the strongly irradiated planets to 
be discussed later, the winds of Jupiter, Saturn, Uranus and Neptune, are significantly slower than 
the planet itself under its own spin (from 12.2 km s^^ for Jupiter to 2.6kms~^ for Neptune, at the 
equator). 

The observed surface winds are believed to be related to motions in the planets' interiors, which, 
according to the Taylor-Proudman theore m, should be confined by the rapid rotation to the plane 
perpendicular to the axis of rotation {e.g. Busse . 19781 ). Unfortunately, no convincing model is yet 
capable of modeling with sufficient accuracy both the interior and the surface layers. 

Our giant planets also exhibit planetary-scale to small-scale storms with very different temporal 
variations. For example, Jupiter's great red spot is a 1 2 000 k m-diameter anticyclone found to have 
lasted for at least 300 years (e.g. Simon-Miller et all 20021). Sto r ms d eveloping over the entire 



planet have even been observed on Saturn ("Sanchez-L avega et al. . 1996f) . Uranus a nd Neptune's 



storm system has been sho wn to have been significantly altered since the Voyager era ([Rages et al 
2OO2I : iHammel et al.l . l2005l). On Jupiter small-scale storms relate d to cumulus-type cloud systems 



have been observed (e.g. Gierasch et all 2000l : iHueso et al. . 2002[ ). and lightning strikes have been 
monitored by Galileo (e.g. lLittle et al.l . ll999l ). These represent only a small arbitrary subset of the 
work concerning the complex atmospheres of these planets. 

It is tempting to extrapolate these observations to the objects outside our Solar System as well. 
However, it is important to stress that an important component of the variability in the atmospheres 
of our giant planets is the presence of relatively abundant condensing chemical species: ammonia 
and water in the case of Jupiter and Saturn, and methane for Uranus and Neptune. These species 
can only condense in very cold atmospheres, thus providing latent heat to fuel important storms. 
Depending on their temperatures and composi t ions, e xtrasolar planets may or may not possess 
such important condensing species re.g. lGuillotl . ll999bh . 



2.7 Moons and rings 

A discussion of our giant planets motivated by the opportunity to extrapolate the results to objects 
outside our solar system would be incomplete without mentioning the moons and rings that these 
planets all possess (see chapters by Breuer & Moore, by Peale and by Husmann et al.). First, the 
satellites/moons can be distinguished from their orbital characteristics as regular or irregular. The 
first ones have generally circular, prograde orbits. The latter tend to have eccentric, extended, 
and/or retrograde orbits. 

These satellites are numerous: After the Voyager era, Jupiter was known to possess 16 satellites, 
Saturn to have 18, Uranus 20 and Neptune 8. Recent extensive observation programs have seen 
the number of satellites increase considerably, with a growing list of satellites presently reaching 
62, 56, 27 and 13 for Jupiter, Saturn, Uranus and Neptune, respectively. All of the new satellites 
discovered since Voyager are classified as irregular. 

The presence of regular and irregular satellites is due in part to the history of planet formation. 
It is believed that the regular satellites have mostly been formed in the protoplanetary subncbulac 
that surrounded the giant planets (at least Jupiter and Saturn) at the time when they accreted 
their envelopes. On the other hand, the irregular satellites are thought to have been captured by 
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the planet. This is for example believed to be the case of Neptune's largest moon, Triton, which 
has a retrograde orbit. 

A few satellites stand out by having relatively large masses: it is the case of Jupiter's lo, 
Europa, Ganymede and Callisto, of Saturn's Titan, and of Neptune's Triton. Ganymede is the 
most massive of them, being about twice the mass of our Moon. However, compared to the mass 
of the central planet, these moons and satellites have very small weights: 10"'^ and less for Jupiter, 
1/4000 for Saturn, 1/25000 for Uranus and 1/4500 for Neptune. All these satellites orbit relatively 
closely to their planets. The farthest one, Callisto revolves around Jupiter in about 16 Earth days. 

The four giant planets also have rings, whose material is probably constantly rcsupplicd from 
their satellites. The ring of Saturn stands out as the only one directly visible with binoculars. 
In this particular case, its enormous area allows it to reflect a sizable fraction of the stellar flux 
arriving at Saturn, and makes this particular ring as bright as the planet itself. The occurrence of 
such rings would make the detection of extrasolar planets slightly easier, but it is yet unclear how 
frequent they can be, and how close to the stars rings can survive both the increased radiation and 
tidal forces. 

2.8 Extrasolar planets 
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Figure 4: Masses and orbital distances of the extrasolar planets discovered by 2006. The size of 
the symbols is proportional to the mass of the parent star (from 0.1 to 4 stellar masses). The color 
(from white to red and black) is proportional to the stellar metallicity. Stars with metallicities 
[Fe/H] < 0.2 are shown in white. The radial velocimetry thresholds at 1 and 10 m/s, encompassing 
the detection limits of most current surveys, are indicated as dashed lines. Transiting planets are 
highlighted with white crosses. 



13 



Huge progresses have been made in the fi eld of extrasolar planets s ince the detection of the 
first giant planet orbiting a solar-type star by iMavor and Quelod (|l995[) . More than 200 planets 
are known at the time of this review, and importantly, 14 planets that transit their star at each 
orbital revolution have been identified (see figS]). These transiting planets are especially interesting 
because of the possibility to measure both their mass and size and thus obtain constraints on their 
global composition. 

In spite of their particular location just a few stellar radii away from their stars, the transiting 
planets that have been discovered bear some resemblance with thei r Solar System cousins in the 
sense that they are also m ostly made of hydrogen and helium (e.g. Burrows et al. . 2000; Guillotl . 
120051: iBaraffe et all l2005l) . They are, however, much hotter due to the intense irradiation that 
they receive. 

Although obtaining direct informations on these planets represent a great observational chal- 
lenge, several key steps ha ve been accomplished: Atomic Sodium, predicted to be detectable 
( Seager and Sas selov. zOOd) , has indeed been detected by transit spectroscopy around one planet 
(jCharbonneau et al.. . 2002) . Hydro dynamically escaping hydrogen, oxygen and carbon have also 
been detected from the same planet ( Vidal-Madiar et al. . 2003L 2004). The measurement of the sec- 
ondary eclipse of several planets by the Spitzer sa tellite allowed a constraint on spectral and hence 



thermal properties of the planetary atmospheres (jCharbonneau et al. ■ I2OO5I: iDeming et all . [2OO5I ) 



Recently, the light curve of a non-transiting planet was detected in the infrared, also with Spitzer 



provi ding preliminary indications of a strong day /nig ht temperature varia tion (jHarrington et al 
l2006f ) perhaps even larger than predicted (see Showman and GuilloU 2002). 



Obviously, there is a big potential for growth in this young field, and the comparison between 
fine observations made for giant planets in our Solar System and the more crude, but also more 
statistically significant data obtained for planets around other stars promises to be extremely 
fruitful to better understand these objects. 



3 The calculation of interior and evolution models 
3.1 Basic equations 

The structure and evolution of a giant planet is governed by the following hydrostatic, thermody- 
namic, mass conservation and energy conservation equations: 



dP 

dr 
dT 
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dL 
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where P is the pressure, p the density, and g ~ Gralr"^ the gravity (m is the mass, r the radius 
and G the gravitational constant). The temperature gradient Vt = (dlnT/dlnP) depends on the 
process by which the internal heat is transported. L is the intrinsic luminosity, t the time, S the 
specific entropy (per unit mass), and e accounts for the sources of energy due e.g. to radioactivity 
or more importantly nuclear reactions. Generally it is a good approximation to assume e for 
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objects less massive than ~ 13 Mj, i.e. too cold to even bum deuterium (but we will see that in 
certain conditions this term may be useful, even for low mass planets). 

The boundary condition at the center is trivial: r = 0; (to — 0, L ^ 0). The external boundary 
condition is more difficult to obtain because it depends on how energy is transpor ted in the atmo- 
spher e. One possibility is to use the Eddington approximation, and to write (e.g. IChandrasekhar , 
I939I ): r — R; {To = TeS, Pq = 2/3 g/n), where Tcff is the effective temperature (defined by 



L = AirRaT^ff, with a being the Stephan-Boltzmann constant), and k is the opacity. Note for 
example that in the case of Jupiter Toff — 124 K, g — 26ms~^ and k « 5 x 10~'^(P/1 bar) m^ kg~^. 
This imphes Pq ~ 0.2 bar (20,000 Pa), which is actually close to Jupiter's tropopause, where 
T « llOK. 

More generally, one has to use an atmospheric model relating the temperature and pressure at a 
given level to the radius R, intrinsic luminosity L and incoming stellar luminosity i*p: r = R; (Tq — 
To{R, L, L^p), Pq = Po{R, L, L^p)). Po is chosen to satisfy the condition that the corresponding 
optical depth at that level should be much larger than unity. If the stellar flux is absorbed mostly 
in a convective zone, then the problem can be simplified by using To{R, L, L^p) fti Tq(R , L + L*p, 0) 
(e.g. Hubbardl. 19771) . An example of such a model is described bv ISaumon et aL (|l996l ) and 



Hubbard et al.l (|2002f ) and is used hereafter to model the planets in the low irradiation limit. 



3.2 High pressure physics & equations of state 

In terms of pressures and temperatures, the interiors of giant planets lie in a region for which accu- 
rate equations of state (EOS) are extremely difficult to calculate. This is because both molecules, 
atoms, and ions can all coexist, in a fluid that is partially degenerate (free electrons have energies 
that are determined both by quantum and thermal effects) and partially coupled (coulomb inter- 
actions between ions are not dominant but must be taken into account). The presence of many 
elements and their possible interactions further complicate matters. For lack of space, this section 
will mostly focus on hydrogen whose EOS has seen the most important developments in recent 
years. A phase diagram of hydrogen (fig. [5]) illustrates some of the important phenomena that 
occur in giant planets. 

The photospheres of giant planets are generally relatively cold (50 to 3000 K) and at low pressure 
(0.1 to 10 bar, or lO'^ to 10^ Pa), so that hydro gen is in molecular form and the perfect gas conditions 
apply. As one goes deeper into the interior hydrogen and helium progressively become fluid. (The 
perfect gas relation tends to underestimate the pressure by 10% or more when the density becomes 
larger than about 0.02 gcm^"^ (P > 1 kbar in the case of Jupiter)). 

Characteristic interior pressures are considerably larger however: as implied by Eqs. [4] and [H 
Pc « GM^/R*, of the order of 10-100 Mbar for Jupiter and Saturn. At these pressures and the 
corresponding densities, the Fermi temperature Tp is larger than 10^ K. This implies that electrons 
are degenerate. Figure [5] shows that inside Jupiter, Saturn, the extrasolar planet HD 209458 b, but 
also for giant planets in general for most of their history, the degeneracy parameter 9 — T/T-p 
is between 0.1 and 0.03. Therefore, the energy of electrons in the interior is expected to be only 
slightly larger than their non-relativistic, fully degenerate limit: Uc > 3/5 fcTp — 15.6 {p/ ^cf^^ eV, 
where k is Boltzmann's constant, /Xc is the number of electrons per nucleon and p is the density in 
gcm~^. For pure hydrogen, when the density reaches ~ O.Sgcm"'^, the average energy of electrons 
becomes larger than hydrogen's ionization potential, even at zero temperature: hydrogen pressure- 
ionizes and becomes metallic. This molecular to metallic transition occurs near Mbar pressures, but 
exactly how this happens remains unclear because of the complex interplay of thermal, coulomb, 
and degeneracy effects (in particular, whether hydrogen metallizes into an atomic state — as 
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Figure 5: Phase diagram for hydrogen with the main phase transitions occurring in the fluid 
or gas phase. The temperature-pressure profiles for Jupiter, Saturn, Uranus, Neptune, and the 
exoplanet HD 209458 b are shown. The dashed nearly vertical line near 1 Mbar is indicative of 
the molecula r to metallic transiti on (here it represents the so-called plasma phase transition as 
calcu l ated by Saumon et aP ( 1995f l). The region in which hydrogen is in solid phase ( Datchi et al.L 
20001 : iGregorvanz et all . \20o£ is represented as a ha tched area. The three phases (I, II, III) of 
solid hydrogen are shown fsee lMao and Hemlevl . Il994l ). Values of the degeneracy parameter 9 are 
indicated as dotted lines to the upper right corner of the figure. 
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suggested in Fig.[S] — or first metallizes in the molecular state H2 remains to be clarified). 

Recent laboratory measurements on fluid deuterium have been able to reach pressures above 
> 1 Mbar, and provide new data in a region where the EOS remai ns most uncertain. Gas- 



guns experiments have been able to measure the reshock temperature (jHolmes et al.l . Il995l ). near 
T ^ 5000 K, P ^ 0.8 Mbar, and a rise in the conductivity of mol ecular hydrogen up to T ~ 3000 K, 
P ~ 1.4 Mbar, a sign that metallicity may have been reached ( Weir et al.l . Il996r). The fol l owing 



few years have see n the development of laser-induc ed shock compression (|da Silva et al . 1997 



Collins et al. . 1998), pul sed-power shock compr ession (Knud son et al. . 2004 ). and convergent shock 



wave experiments (Belov et al. . 20021: iBoriskov et al. . .2005,) in a high-pressure (P = 0.3 — 4 Mbar) 
high-temperature (T ~ 6000 — 10^ K) regime. Unfortunately, experimental results along the prin- 
cipal Hugoniot of deuterium do not agree in this pressure range. Laser compression data give 
a maximum compression of ~ 6 while both the pulsed-power compression experiments and the 
convergent shock wave experiments find a value of ^ 4 . Models that ar e partly calibrated with 
experimental data ( Saumon et al. . 1995 : Rosd . 1998; Ross and Yane . 200l) obtain a generally good 
agreement with the laser-compression data. However, the fact that indep endant models based on 
first principles ( Militzer et al.l . 2001; Desiarlai^ , 2003t Bonev et al. . 20041 ) yield low compressions 
strongly favors this solution. 

The question of the existence of a first-order molecular to metallic transition of hydrogen 
(i.e. both molecular dissociation and ionisation occur simultaneously and discontinuously at the 
so-called plasma phase transit ion, or PPT) r emain s however. The critical line shown in fig. [5] 
corresponds to calculations by ISaumon et all but may be caused by artefacts in the free 

energy calculation. Recent Density Functional Theory (DFT) simulations bv iBonev et al. ( 2004h 
i ndicate the possibili ty of a first order liquid-liquid transition but other path-integral calculations 
( Militzer et al.l . [200ll ) do not. It is crucial to assess the existence of such a PPT because it would 
affect both convection and chemical composition in the giant pl anets. 

A clear result from fig. [5] at least is that, as first shown by HubbardI ( 19681 ). the interiors of 
the hydrogen-helium giant planets are fluids whatever their age: of course, they avoid the critical 
point for the liquid gas transition in hydrogen and helium, at very low temperatures, but they also 
lie comfortably above the solidification lines for hydrogen and helium. (An isolated Jupiter should 
begin partial solidification only after at least ^ 10'^ Ga of evolution.) They are considered to be 
fluid because at the high pressures and relatively modest temperatures in their interiors, coulomb 
interactions between ions play an important role in the EOS and yield a behavior that is more 
reminiscent of that of a liquid than that of a gas, contrary to what is the case in e.g. solar-like stars. 
For Uranus and Neptune, the situation is actually more complex because at large pressures they 
are not expected to contain hyd rogen, b ut numerical simulations show that ices in their interior 
should be fluid as well (Cavazzoni et al. I, [1999) ■ 

Models of the interiors of giant planets require thermodynamically consistent EOSs calculated 
over the entire domain of pressure and temperature spanned by the planets during their evolution. 
Elements other than hydrogen, most importantly helium, should be consistently included. Such 
a calculation is a daunting t ask, and the only rec ent attempt at such an astrophysical EOS for 
substellar objects is that by Saumon et al.l ( 1995f ) . Another set of EOSs reproducing either the 
high- or low-compression results was calculated by ISaumon and Gufllot (,2004) specifically for the 
calculation of present-day models of Jupiter and Saturn. 

These EOSs have so far included other elements (including helium), only in a very approximative 
way, i.e. with EOSs for helium and heavy elements that are based on interpolations between 
somewhat ideal regimes, using an additive volume law, and neglect ing the possibility of existence 
of phase separations (see Hubbard et al.l . [2OO2I : iGuillot et al. . 2004L for further discussions). 
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3.3 Heat transport 

Giant planets possess hot interiors, implying that a relatively large amount of energy has to be 
transported from the deep regions of the planets to their surface. This can either be done by radi- 
ation, conduction, or, if these processes are not sufficient, by convection. Convection is generally 
ensured by the rapid rise of the opacity with increasing pressure and temperature. At pressures of 
a bar or more and relatively low temperatures (less than 1000 K), the three dominant sources of 
opacities are water, methane and collision-induced absorption by hydrogen molecules. 

However, in the intermediate temperature range between ^ 1200 and 1500 K, the Rosseland 
opacity due to the hydrogen and helium absorption behaves differently: the absorption at any given 
wavelength increases with density, but because the temperature also rises, the photons are emitted 
at shorter wavelengths, where the monochromatic absorption is smaller. As a consequence, the 
opacity can decrease. This was shown by Guillot et al. (1994) to potentially lead to the presence 
of a deep radiative zone in the interiors of Jupiter, Saturn and Uranus. 

This problem must however be reanalyzed in the light of recent observations and analyses of 
brown dwarfs. Their spectra show unexpectedly wide sodium and potassium absorption lines (see 
Burrows, Marley & Sharp 2000), in spectral regions where hydrogen, helium, water, methane and 
ammonia are relatively transparent. It thus appears that the added contribution of these elements 



(if th ey are indeed present) would wipe out any radiative region at these levels (|Guillot et al 
20041 ). 



At temperatures above 1500 ^ 2000 K two important sources of opacity appear: (i) the rising 
number of electrons greatly enhances the absorption of and ; (ii) TiO, a very strong absorber 
at visible wavelengths is freed by the vaporization of CaTiOs. Again, the opacity rises rapidly which 
ensures a convective transport of the heat. Still deeper, conduction by free electrons becomes more 
efficient, but the densities are foun d not to b e high enou gh for th is process t o be sig nificant, except 



erncient, out the densities are loun g not to b e hign enou gh tor th is process t o be sig i 
perhaps near the central core (see Hubbardl . [l968 : Stevenson and Salpeter , 1977lJ l. 



While our giant planets seem to possess globally convective interiors, strongly irradiated ex- 
trasolar planets must develop a radiative zone just beneath the levels where most of the stellar 
irradiation is absorbed. Depending on the irradiation and characteristics of the planet, this zone 
may extend down to kbar levels, the deeper levels being convective. In this case, a careful deter- 
mination of the opacities is necessary (but generally not po ssible) as these control the cooling and 
contraction of the deeper interior (see lFerguson et aL . 1200,4 for a discussion of opacities and tables 
for substellar atmospheres and interiors). 

3.4 The contraction and cooling histories of giant planets 

The interiors of giant planets is expected to evolve with time from a high entropy, high 9 value, 
hot initial state to a low entropy, low 0, cold degenerate state. The essential physics behind can 
be derived from the well-known virial theorem and the energy conservation which link the planet's 
internal energy iJj, gravitational energy Eg and luminosity through: 

^E, + Eg = 0, (8) 

where f — J^^ 3{P/p)dm/ udm ~<iP/ pu>, the brackets indicating averaging, and u is the 
specific internal energy. For a diatomic perfect gas, ^ = 3.2; for fully-degenerate non-relativistic 
electrons, ^ = 2. 
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Thus, for a giant planet or brown dwarf beginning its life mostly as a perfect H2 gas, two third 
of the energy gained by contraction is radiated away, one third being used to increase Ei. The 
internal energy being proportional to the temperature, the effect is to heat up the planet. This 
represents the slightly counter-intuitive but well kno wn effect that a star or giant pla net initially 



heats up while radiating a significant luminosity (e.g. iKippenhahn and Weigertl . 119941 ). 

Let us now move further in the evolution, when the contraction has proceeded to a point where 
the electrons have become degenerate. For simplicity, we will ignore coulomb interactions and 
exchange terms, and assume that the internal energy can be written as = Ed + Eiom and 
that furthermore E^i ^ -Eion is small). Because ^ « 2, we know that half of the gravitational 
potential energy is radiated away and half of it goes into internal energy. The problem is to decide 
how this energy is split into an electronic and an ionic part. The gravitational energy changes 
with some average value of the interior density as Eg oc 1/R oc p^/^. The energy of the degenerate 
electrons is essentially the Fermi energy: Ed oc p^/^. Therefore, Eci ~ 2(i?ci/£'g)i?g. Using the 
virial theorem, this yields: 

^1 ~ -Eg w 2L (10) 

L « -Eion^-t. (11) 

The gravitational energy lost is entirely absorbed by the degenerate electrons, and the observed 
luminosity is due to the thermal cooling of the ions. 

Several simplifications limit the applicability of this result (that would be valid in the white 
dwarf regime). In particular, the coulomb and exchange terms in the EOS introduce negative 
contributions that cannot be neglected. However, the approach is useful to grasp how the evolution 
proceeds: in its very early stages, the planet is very compressible. It follows a standard Kelvin- 
Helmoltz contraction. When degeneracy sets in, the compressibility becomes much smaller {aT ^ 
0.1, where a is the coefficient of thermal expansion), and the planet gets its luminosity mostly 
from the thermal cooling of the ions. The luminosity can be written in terms of a modified Kelvin- 
Helmoltz formula: 

L.,^, (12) 

where r is the age, and 77 is a factor that hides most of the complex physics. In the approximation 
that coulomb and exchange terms can be neglected, rj m 9/ {9+1). The poor compressibility of giant 
planets in their mature evolution stages imply that 77 <SC 1 (77 0.03 for Jupiter): the luminosity 
is not obtained from the entire gravitational potential, but from the much more limited reservoir 
constituted by the thermal internal energy. Equation 1121 shows that to first order, logL oc — logr: 
very little time is spent at high luminosity values. In other words, the problem is (in most cases) 
weakly sensitive to initial conditions. However, it is to be noticed that with progresses in our 
capabilities to detect very young objects, i.e. planets and brown dwa rfs of only a few mi llion years 



of age, the problem of the initial conditions does become important (jMarlev et al.l . 120061 . 

Figure [S] shows more generally how giant planets, but also brown dwarfs and small stars see 
their luminosities evolve as a function of time. The l/r slope is globally conserved, with some 
variations for brown dwarfs during the transient epoch of deuterium burning, and of course for 
stars, when they begin burning efficiently their hydrogen and settle on the main sequence: in 
that case, the tendency of the star to contract under the action of gravity is exactly balanced by 
thermonuclear hydrogen fusion. 
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Figure 6: Evolution of the luminosity (in Lq) of solar- metallicity M dwarfs and substellar objects 
vs. time (in yr) after formation. In this figure, "brown dwarfs" are arbitrarily designated as those 
objects that burn deuterium, while those that do not are tentatively labelled "planets". Stars are 
objects massive enough to halt their contraction due to hydrogen fusion. Each curve is labelled by 
its corresponding mass in Mq, with the lowest t hree correspond i ng to the mass of Saturn, half the 
mass of Jupiter, and the mass of Jupiter. [From iBurrows et al.l (|l997l )]. 
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Figure 7: Theoretical and observed mass-radius relations. The black line is applicable to the 
evolution of solar composition planets, brown dwarfs and stars, when isolated or nearly isolated 
(as Jupiter, Saturn, Uranus and Neptune, defined by diamonds and their respective symbols), after 
5 Ga of evolution. The dotted line shows the effect of a 15 core on the mass-radius relation. 
Orange and yellow curves represent the mass-radius relations for heavily irradiated planets with 
equilibrium temperatures of 1000 and 2000 K, respectively, and assuming that 0.5% of the incoming 
stellar luminosity is dissipated at the center (see section 14. 3p . For each irradiation level, two cases 
are considered: a solar-composition planet with no core (top curve), and one with a 100 M0 
central core (bottom curve) . The transiting extrasolar giant planets for which a mass and a radius 
was measured are shown with points that are color-coded in function of the planet's equilibrium 
temperature. The masses and radii of very low mass stars are also indicated as blue points with 
error bars. 
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3.5 Mass-radius relation 



The relation between mass and radius has very fundamental astrophysical applications. Most 
importantly it allows one to infer the gross composition of an object from a measurement of its 
mass and radius. This is especially relevant in the context of the discovery of extrasolar planets 
with both radial velocimetry and the transit method, as the two techniques yield relatively accurate 
determination of M and R. 

Figure [7] shows mass-radius relations for compact degenerate objects from giant planets to 
brown dwarfs and low- mass stars. The right-hand side of the diagram shows a rapid increase of 
the radius with mass in the stellar regime which is directly due to the onset of stab le thermonuc lear 
reactions. In this regime, observations and theoretical models agree (see however Ribas . 20061 for 
a more detailed discussion). The left-hand side of the diagram is obviously more complex, and 
this can be understood by the fact that planets have much larger variations in compositions than 
stars, and because external factors such as the amount of irradiation they receive do affect their 
contraction in a significant manner. 

Let us first concentrate on isolated or nearly-isolated gaseous planets. The black curves have a 
local maximum near 4Mj: at small masses, the compression is small so that the radius increases 
with mass. At large masses, degeneracy sets in and the radius decreases with mass. 

This can be understood on the basis of polytropic models based on the assumption that P = 
^pi+i/"^ where K and n are constants. Because of degeneracy, a planet of large mass will tend 
to have n 1.5, while a planet a smaller mass will be less compressible {n — > 0). Indeed, it can 
be shown that in their inner 70 to 80% in radius isolated solar composition pla nets of 10, 1 and 
0.1 M .T have n = 1.3, 1.0 and 0.6, respectively. From polytropic equations (e.g. IChandrasekliar , 
19391) : 

RckK^M^. (13) 

Assuming that K is independant of mass, one gets R oc M^ -*^^, Af°, and M~°-^^ for M = 10, 1 
and 0.1 Mj, respectively, in relatively good agreement with fig. [7] (the small discrepancies are due 
to the fact that the intrinsic luminosity and hence K depend on the mass considered). 

Figure [7] shows already that the planets in our Solar System are not made of pure hydrogen 
and helium and require an additional fraction of heavy elements in their interior, either in the form 
of a core, or distributed in the envelope (dotted hue). 

For extrasolar planets, the situation is complicated by the fact that the intense irradiation that 
they receive plays a major role in their evolution. The present sample is already quite diverse, 
with equilibrium temperature (defined as the effective temperature corresponding to the stellar flux 
received by the planet) ranging from 1000 to 2500 K. Their composition is also quite variable , with 
some planets having large masses of heavy elements (jSato et al.l . 1200,4 iGuillot et al.l . l2006h . The 
orange and yellow curves in fig.[7]show theoretical results for equilibrium temperatures of 1000 and 
2000 K, respectively. Two extreme models have been plotted: assuming a purely solar composition 
planet (top curve), and assuming the presence of a 100 central core (bottom curve). In each 
case, an additional energy source proportional to 0.5% of the incoming luminosity was also assumed 
(see discussion in § 14.31 hereafter). 

The increase in radius for decreasing planetary mass for irradiated, solar-composition planets 
with little or no core can be understood using the polytropic relation (eq. I13p . but accounting 
for variations of K as defined by the atmospheric boundary condition. Using the Eddington 
approximation, assuming k oc P and a perfect gas relation in the atmosphere, one can show 
that K oc (M/i?2)-i/2« and that therefore R oc M^^. With n = 1, one finds R oc M-^/'^. 
Strongly irradiated hydrogen-helium planets of small masses are hence expected to have the largest 
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radii which qualitatively explain the positions of the extrasolar planets in fig. [71 Note that this 
estimate implicitly assumes that n is constant throughout the planet. The real situation is more 
complex because of the growth of a deep radiative region in most irradiated planets, and because 
of structural changes between the degenerate interior and the perfect gas atmosphere. 

In the case of the presence of a fixed mass of heavy elements, the trend is inverse because of the 
increase of mean molecular mass (or equivalently core/envelope mass) with decreasing total mass. 
Thus, small planets with a core are much more tightly bound and less subject to evaporation than 
those that have no core. 



3.6 Rotation and the figures of planets 

The mass and radius of a planet informs us on its global composition. Because planets are also 
rotating, one is allowed to obtain more information on their deep interior structure. The hydrostatic 
equation becomes more complex however: 

^^fr^ fff P(0 ,3 



p ^^{^111 |„ ] - ft X {fl X r), (14) 



where ft is the rotation vector. The resolution of eq. (|14p is a complex problem. It can however 
be somewhat simplified by assuming that |r2| = a; is such that the centrifugal force can be derived 
from a potential. The hydrostatic equilibrium then writes VP = pVt/, and the figure of the 
rotating planet is the n defined by the U = constant level surface. 

One can show (e.g. Zharkov and Trubitsvn . 19781 ) that the hydrostatic equation of a fluid planet 



can then be written in terms of the mean radius r (the radius of a sphere containing the same 
volume as that enclosed by the considered equipotential surface): 

IdP Gm 2 2_ GM_ 

= — =2- + ^ + ^3-'''/''^' (15 
p or r 3 

where M and R are the total mass and mean radius of the planet, and fuj is a slowly varying 
function of r. (In the case of Jupiter, ip^ varies from about 2 x 10"'^ at the center to 4 x 10^"^ at 
the surface.) Equations (IS][7]) remain the same with the hypothesis that the level surfaces for the 
pressure, temperature, and luminosity are equipotentials. The significance of rotation is measured 
by the ratio of the centrifugal acceleration to the gravity: 

As discussed in section [2?2l in some cases, the external gravity field of a planet can be accu- 
rately measured in the form of gravitational moments Jk (with zero odd moments for a planet in 
hydrostatic equilibrium) that measure the departure from spherica l symmetry. Together w i th the 
mass, this provides a constraint on the interior density profile (see Zharkov and Trubitsvnl ( 1974[ ) 
-see also chapters by Van Hoolst and Sohl & Schubert): 



M - JJJ pir,e)d-'T, 

^2, = JJJ Pir,ey^P2.dcos0)d\, 

where dr is a volume element and the integrals are performed over the entire volume of the planet. 
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Figure [5] shows how the different layers inside a planet contribute to the mass and the grav- 
itational moments. The figure applies to Jupiter, but would remain relatively similar for other 
planets. Note however that in the case of Uranus and Neptune, the core is a sizable fraction of 
the total planet and contributes both to J2 and J4. Measured gravitational moments thus provide 
information on the external levels of a planet. It is only indirectly, through the constraints on the 
outer envelope that the presence of a central core can be infered. As a consequence, it is impossible 
to determine this core's state (liquid or solid), structure (differentiated, partially mixed with the 
envelope) and composition (rock, ice, heliimr...). 




R/Rtot 



Figure 8: Contribution of the level radii to the gravitational moments of Jupiter. Jq is equivalent 
to the planet's mass. The small discontinuities are caused by the following transitions, from left 
to right: core/envelope, helium rich/helium poor (metallic/molecular). Diamonds indicate the 
median radius for each moment. 



For planets outside the solar system, although measuring their gravitati onal potential is utopic , 
their oblateness may be reachable with future space transit observations (jSeager and Huil . 120021 ). 
Since the oblateness e is, to first order, proportionnal to q: 



R. 



cq 



Rpol 



1 



(17) 



(where A2 = J2/Q ~ 0.1 to 0.2), it may be possible to obtain their rotation rate, or with a rotation 
measured from another method, a first constraint on their interior structure. 



4 Interior structures and evolutions 
4.1 Jupiter and Saturn 

As illustrated by fig.El the simplest interior models of Jupiter and Saturn matching all observational 
constraints assume the presence of three main layers: (i) an outer hydrogen-helium envelope, whose 
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Ices + Rocks 
core ? 

Jupiter Saturn 

Figure 9: Schematic representation of the interiors of Jupiter and Saturn. The range of tempera- 
tures is estimated using homogeneous models and including a possible radiative zone indicated by 
the hashed regions. Helium mass mixing ratios Y are indicated. The size of the central rock and ice 
cores of Jupiter and Saturn is very uncertain (see text). In the case of Saturn, the inhomogeneous 
region may extend down all the way to the core which would imply the formation of a helium core. 
[Adapted fromiGuUlot (,1999b^ ]. 
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global composition is that of the deep atmosphere; (ii) an inner hydrogen-helium envelope, enriched 
in helium because the whole planet has to fit the H/He protosolar value; (iii) a central dense 
core. Because the planets are believed to be mostly convective, these regions are expected to be 
globally homogeneous. (Many interesting thermochemical transformations take place in the deep 
atmosphere, but they are of little concern to us). 

The transition from a helium-poor upper envelope to a helium-rich lower envelope is thought 
to take place through the formation of helium-rich droplets that fall deeper into the planet due 
to their larger density. These droplets form because of an assumed phase transition of helium 
in hydrogen at high pressures and low temperatures. Three-layer models implicitely make the 
hypothesis that this region is narrow. Indeed, calculations of such a phase separation in a fully- 
ionized plasma indicate a rapid decrease of the critical temperature with increasing pressure, with 
the consequence that helium would be unsoluble in a relatively s mall region i n low- pressure metallic 
hydrogen. This region would progressively grow with time (e.g. Stevenson . 19821 ). However, DFT 
calculation s have indicate d that the critical temperature for helium demixing may rise with pressure 
(jPfaffenzcUer et al.l . [l995l ) . presumably in the regime where hydrogen is only partially ionized and 
bound states remain. This opens up the possibility that the inhomogeneous regions may be more 
extended, and that models more complex than the three-layer models may be needed, in particular 
in the case of Saturn (see below). 

In the absence of these calculations, the three-layer models can be used as a useful guid- 
ance to a necessarily hypothetical ensemble of allowed structures and compositions of Jupiter 
and Saturn. These relative ly extensive exploration of the parameter space have been performed by 
Saumon and GiiillotJ (j2004h . The calculations assume that only helium is inhomogeneous in the en- 



velope (the abundance of heavy elements is supposed to be uniform accross the molecular/metallic 
hydrogen transition). Many sources of uncertainties are taken into account however; among them, 
the most significant are on the equations of state of hydrogen and helium, the uncertain values of J4 
and Jfi, the presence of differential rotation deep inside the planet, the location of the helium-poor 
to helium-rich region, and the uncertain helium to hydrogen protosolar ratio. 

Their results indicate that Jupiter's core is smaller than ^ 10 M^, and that its global com- 
position is pretty much unknown (between 10 to 42 M^ of heavy elements in total). The models 
indicate that Jupiter is enriched compared to the solar value by a factor 1.5 to 8 times the solar 
value. This enrichment is compatible with a global uniform enrichment of all species near the 
atmospheric Galileo values, but include many other possibilities. 

In the case of Saturn, the solutions depend less on the hydrogen EOS because the Mbar pressure 
region is comparatively smaller. The total amount of heavy elements present in the planet can 
therefore be estimated with a better accuracy than for Jupiter, and is between 20 and 30 M^. In 
three-layer models with a discontinuity of the helium abundance at the molecular-metallic hydrogen 
interface but continuity of all other elements, the core masses found are between 10 and 22 M^. 
However, because Saturn's metallic region is deeper into the planet, it mimics the effect that a 
central core would have on J2. If we allow for variations in the abundance of heavy elements 
together with the helium discontinuity, then the core mass can become much smaller, and even 
solutions with no core can be found (Guillot 1999a). These solutions depend on the hypothetic 
phase separation of an abundant species (e.g. water), and generally cause an energy problem 
because of the release of considerable gravitational energy. However, another possibility is through 
the formation of an alm ost pure heliu m shell around the central core, which could lower the core 
masses by up to TMg (IFortnev and H ubbard. 2003, Hubbard, personnal communication). 

Concerning the evolutions of Jupiter and Saturn, the three main sources of uncertainty are, 
by order of importance: (1) the magnitude of the helium separation; (2) the EOS; (3) the at- 



26 



mospheric boundary conditions. Figure [TO] shows an ensemble of possibilities that attempt to 
bracket the minimum and maximum c ooling. In all cases, helium sedimentation is needed to ex- 
plain Saturn's present luminosity (see ISalpeteii . Il973t IStevenson and Salpeten . Il977at iHubbard , 
19771 ). Recent models of Saturn 's evolution appear to favor a scenario in which h elium settles 
down almost to the central core ( Hubbard et al. . 119991 : iFortney and Hubbardl - liooal ). In the case 
of Jupiter, the sedimentation of helium that appears to be necessary to explain the low atmospheric 
helium abundance poses a problem for evolution models because it appears to generally prolong its 
evolution beyond 4.55 Ga, the age of the Solar System. However, different solutions are possible, 
including improvements of the EOS and atmospheric boundary conditions, or even the possible 
progressive erosion o f the central core that would yield a lower Jupiter's luminosity at a given age 



(iGuillot et all . 120041) 
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Figure 10: Final stages of evolution of Jupiter and Saturn. The present effective temperatures, 
reached after ~ 4.55 Ga of evolution, arc indicated as horizontal orange lines. For each planet 
two models represent attempts to bracket the ensemble of possibilities, with the faster evolu- 
tion corresponding to that of an homogeneous planet, while the slo west evolution include s the 
effect of helium settling in the last evolution phase. [Adapted from Hubbard et al. ( 1999[ ) and 
Fortnev and HubbardI (|2003l )]. 



4.2 Uranus and Neptune 

Although the two planets are relatively similar, fig. [7] already shows that Neptune's larger mean 
density compared to Uranus has to be due to a slightly different composition: either more heavy 
elements compared to hydrogen and helium, or a larger rock/ice ratio. The gravitational moments 
impose that the density profiles lie close to that of "ices" (a mixture initially composed of e.g. 
H2O, CH4 and NH3, but which rapidly becomes a ionic fluid of uncertain chemical composition in 
the planetary interi or), except in the outermost layers , which have a density closer to that of hy- 
drogen and helium ( Marlev et al.l . 1995t Podolak et al. . 2000f l. As illustrated in fig. [HI three-layer 
models of Uranus and Neptune consisting of a central "rocks" core (magnesium-silic ate and iron 
material), an ice layer and a hydrogen-helium gas envelope have been calculated (|Podolak et all 
199ll : lHubbard et al.l . lT995l) . 
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Rocks? 



Uranus 



Neptune 



Figure 11: Schematic representation of the interiors of Uranus and Neptune. [Adapted from lGuillotI 
(|l999bh ]. 



The fact that models of Uranus assuming homogeneity of each layer and adiabatic temperature 
profiles fail in reproducing its gravitational mo ments seem to imply that substantial parts of the 
planetary interior are not homogeneously mixed ( Podolak et al. . [19951 ). This could explain the fact 
that Uranus' heat flux is so small: its heat would not be allowed to escape to space by convection, 
but through a much slower diffusive process in the regions of high molecular weight gradient. 
Such regions would also be present in Neptune, but much deeper, thus allowing more heat to be 
transported outward. The existence of these non-homogeneous, partially mixed regions are further 
confirmed by the fact that if hydrogen is supposed to be confined solely to the hydrogen-helium 
envelope, models predict ice/rock ratios of the order of 10 or more, much larger than the protosolar 
value of ~ 2.5. On the other hand, if we impose the constraint that the ice/rock ratio is protosolar, 
the overall composition of both Ura nus and Neptune is, by mass, about 25% rock, 6 — 70% ice, 
and 5 - 15% hydrogen and helium (jPodolak et all . Il99ll Il995t iHubbard et"al] . Il995f) . Assuming 
both ice and rock are present in the envelope, an upper limit to the amount of hyd rogen and helium 



present is ~ 4.2 Me for Uranus and ~ 3.2 for Neptune (jPodolak et all . |2000| ) . A lower limit of 
^ 0.5 Mg for both planets can be inferred by assuming that hydrogen and helium are only present 
in the outer envelope at P < lOOkbar. 



4.3 Irradiated giant planets 

Although all extrasolar giant planets are in principle interesting, we focus here on the ones that orbit 
extremely close to their star because of the possibility to directly characterise them and measure 
their mass, radius and some properties of their atmosphere. Two planets are proxies for this new 
class of objects: the first extrasolar giant planet discovered, 51 Peg b, with an orbital period of 
P = 4.23 days, and the first transiting extrasolar giant planet, HD 209458 b, with P = 3.52 days. 
Both planets belong to the Pegasus constellation, and following astronomical conventions (e.g. 
Cepheids, named after 6 Cephei), we choose to name giant planets orbiting close to their stars 
with periods shorter than 10 days "Pegasids" (alternatively, "hot Jupiters" is also found in the 
litterature) . 

With such a short orbital period, these planets are for most of them subject to an irradiation 
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from their central star that is so intense that the absorbed stellar energy flux can be about ~ 10^ 
times larger than their intrinsic flux. The atmosphere is thus prevented from cooling, with the 
consequence that a r adiative zone develops and governs the cooling and contraction of the interior 



(|Guillot et all . 119961 ). Typically, for a planet like HD 209458 b, this radiati ve zone extends to kbar 



levels , T ^ 4000 K, and is located in the outer 5% in radius (0.3% in mass) ( Guillot and ShowmanL 
20021 ). 




Figure 12: Conjectured dynamical structure of Pegasids (strongly irradiated extrasolar giant plan- 
ets): At pressures larger than 100-800 bar, the intrinsic heat flux must be transported by convec- 
tion. The convective core is at or near synchronous rotation with the star and has small latitudinal 
and longitudinal temperature variations. At lower pressures a radiative envelope is present. The 
top part of the atmosphere is penetrated by the stellar light on the day side. The spatial variation 
in insolation should d r ive w inds that transport heat from the day side to the night side. [From 
Showman and GuillotI (|2002l )]. 



Problems in the modeling of the evolution of Pegasids arise mostly because of the uncertain 
outer boundary condition. The intense stellar flux implies that the atmospheric temperature profile 
is extremely dependant upon the opacity sources considered. Depending on the chosen composi- 
tion, the opacity data used, the assumed presence of clouds, the geo metry considered, re s ulting 



temperatu res in the deep atmosphere can differ by up to ~ 600 K (iSe ager and Sasselo vl 12000 ^ 



Goukenleuaue et al. I. l2000l :l iBarman et al.l . [2001 :I Sudarskv et al.l 1. 120031: 1 liro et al.. , .2005 : Fortnev et al 



20061 ). Furthermore, as illustrated by fig. [TH the strong irradiation and expected synchronisa- 



tion of the planets implies that strong inhomogeneities should exist in the atmosphere with in 
particular strong (~ 500 K) day-night and equator-to-pole differences in effective temperature s 
(jShowman and GuillotI . l2002l : liro et all . 120051 : ICooper and Showmanl . l2005t iBarman et all . |2005| ). 
further complicating the modeling of the planetary evolution (see fig. fig. [T5|) . Finally, another 
related problem is the presence of the radiative zone. Again, the composition is unknown and 
the opacity data are uncertain in this relatively high temperature (T ~ 1500 — 3000 K) and high 
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pressure (up to ^ 1 kbar) regime. 

We have seen in fig. fig. [7] that the measured masses and radii of transiting planets can be 
globally explained in the framework of an evolution model including the strong stellar irradiation 
and the presence of a variable mass of heavy elements, either in the form of a central core, or 
spread in the planet interior. However, when analyzing the situation for each planet, it appears 
that several planets are too large to be reproduced by standard models, i.e. models using the 
most up-to-date equations of state, opacities, atmospheric boundary conditions and assuming that 
the planetary luminosity governing its cooling is taken solely from the lost gravitational potential 
energy (see Section [XTt . 

Figure [HI illustrates the situation for the particular case of HD209458b: unless using an unre- 
alistically hot atmosphere, or arbitrarily increasing the internal opacity, or decreasing the helium 
content, one cannot reproduce the observed radius which i s 10 to 20% larger than calculated 



20031) . The fact that 



(jBodenheimer et al.l . l200ll l2003t iGuillot and Showmanl l2002t iBaraffe et al. 

the measured radius corresponds to a low-pres sure (^mbar) le vel while the calculated radius cor- 
responds to a level near Ibar is not negligible ( Burrows et al.l . [2003 ) but too small to account for 
the difference. This is problematic because while it is easy to invoke the presence of a massive core 
to explain the small size of a planet, a large size such as that of HD209458b requires an additional 
en ergy source, or signif i cant m odifications in the data/physics involved. 

Bodenheimer et all ( 2001 ) proposed that this large radius may be due to a small forced ec- 



centricity (e ^ 0.03) of HD209458b, and subsequent tidal dissipation in the planet inter ior, but 
detai led observations indicate that the eccentricity is small, e — 0.014 ± 0.009 (Laughlin et al 



20051), and observations o f the secondary eclip se imply that this would further require a chance 



configuration of the orbit (jPeming et al. Another proposed explanation also involving tidal 



dissipation of or bital energy is that t he pl anet may be trapped in a Cassini state with a large or- 
bital inclination ( Winn and Holmanl . 1200,^ . but it appears to have a low probability of occurrence 
( Levrard et al. . 2006h . Finallv. a third possibility that would apply to all Pegasids is to invoke a 
downward transport of kinetic energy and its dissipation by tides ( Showman and Guillotl . 2002 ). 
This last possibility would re quire the vario us transiting planets to have different core masses to 
reproduce the observed radii (|Guillotl . [2005h . 

Recently, as more transiting Pegasids have been discovered, the number of anomalously large 
ones has increased to at least 3 for 11 planets, implying that this is not a rare event. This lends 
more weight to a mechanism that would apply to each planet. In this case, masses of heavy 
elements can be derived by imposing that all planets should be fitted by the same n iodel with the 
same hypotheses. This can be done by inverting the results of fig. [3 as described by ( Guillot et al.l 
20061 ). The method is applied to the known transiting Pegasids by the end of 2006 in fig. [151 a 



plot of the masses of heavy elements in the planets as a function of the metallicities of the parent 
star (which measures how rich a given star is in heavy elements compared to the Sun). 

Figure [TSl first shows that in some cases, large masses of heavy elements (up to ^ 100 M0 are 
necessary. This is in harmony with the composition inferred for HD149026b, i.e. around 70 
of heavy elements, a conclus ion that is hard to escape because of t he low total mass and high 
irradiation of the planet (see llkoma et al. ■ l2006t iFortnev et al.l . [2006l ). Furthermore, there seems 
to be a correlation between th e mass of heavy elem ents inferred in the transiting planets, and the 
metallicity of the parent stars ( Guillot et al. . 20061 ). although this correlation has to be ascertained 
by more measurements. A caveat is important: these results are intrinsically model-dependent 
as they are based on the assumption that all planets receive an additionnal tidal heat flux that 
is proportionnal to the energy flux that they receive in the form of photons from the star. More 
transiting planets are needed to confirm or infirm this model, but the large variety in core masses, 
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Figure 13: Temperature versus pressure for a sequence of locations in the atmosphere of 
HD209458b, assuming no horizontal redistribution of heat. Each sequence corresponds to a given 
direction of the incident flux relative to the surface normal. The approximate regions represented 
by the collection of T-P profiles are shown as solid black lines on the illustrative sphere. The top- 
most T-P profile corresponds to the sub-stellar point (black dot on the sphere). The terminator 
and night side (black hemisphere) are modeled with the non- irradiated profile (lowest T-P curve). 
The radiative-convective boundary at the sub-stellar point and on the night side are labeled with 
filled circles. The dashed lines indicate the approximate condensation curves for three common 
grain species. The dotted line indicates where gaseous CO and CH4 concentrations are equal (CO 
is dominant to the left of this line). The thick, grey, dashed lines are T-P profiles calculated for 
a normal incident fiux equal to 0.5 (top) and 0.25 (bottom) times that at the substellar point, as 
often used as approx imate solutions for the day side, or entire atmosphere, respectively. [From 



Barman et all (|2005l )]. 
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Figure 14: The contraction of HD209458b as a function of time can be compared to its measured 
radius and inferred age shown by the black box. Standard models (blue curve) for the evolution of 
that 0.69 Mj planet generally yield a radius that is too small compared to the observations, even 
for a solar composition and no central core (a larger core and -in most cases- larger amounts of 
heavy elements in the planet imply an even smaller size for a given age) . Unrealistically low helium 
abundances or high opacities models lead to evolution tracks that barely cross the observational 
box. A possiblity is that heat is dissipated into the deep interior by stellar tides, either related to 
a non-zero orbital eccentricity forced by an unseen companion, or because of a constant transfer 
of angular momentum from the heated atmosphere to the interior (black curve). Alternatively, 
the atmosphere may be hotter than predicted due to heating by strong zonal winds and shear 
instabilities (red curve). 
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Figure 15: Mass of heavy elements in transiting Pegasids known by 2006 as a function of the metal 
content of the parent star relative to the Sun. The mass of heavy elements required to fit the 
measured radii is calculated on the basis of evolution models including an additional heat source 
slowing the cooling of the planet. This heat source is assumed equal to 0.5% of the incoming stellar 
heat flux fShowman and Guillotl . l20ol ^. Horizontal error bars correspond to the la errors on the 
[Fe/H] determination. Vertical error bars are a consequence of the uncert a inties on the measured 
planetary radii and ages. Note that the results, based on Guillot et al. ( 20061 ) are intrinsically 
model-dependent and may be affected by further discoveries of transiting planets. 
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and the absence of Pegasids around metal-poor stars ([Fe/H] < —0.07) as indicated by fig. [T5] 
appear to be robust consequences of this work. 

Another intriguing possibility concerning Pegasids is that of a sustained mass loss due to the 
high i rradiation dose that the planets receive. I ndeed, this effe ct was predicted ( Burrows and Lunin3 . 
1995HGuillot et aLLlToollLammer eraD . l2003l ) and detected (fVidal-Madjar et all. 1 2003'. '2004 ), but 
its magnitude is still quite uncertain, by at le ast two orders of magnitude (jLammer et al. . 2003 : 



Lecavelier des Etangs et al. . 2004; Yellel . 20061) . The effect on the evolution is surprisingly lim- 



ited, except a. t the final stages when an exponential mass loss appears in fully gaseous planets 



( Baraffe et al.U2004 ) 



Finally, it is important to note that another class of planets awaits a direct characterisation by 
I transit method: that of ice or rock giants. Small-mass planets around 10 have been detected 





Lovis et al. 


,2006 


; Beaulieu et al. 


1996 


: Valencia et al.. 


20061), and we c 



2ooir y not have the observational capability to test 

whether they transit in front of their star. This should be resolved by the space mission CoRoT 
(launched on 27 dec 2006) and Kepler (launch ^ 2008). These objects are especially interesting 
but pose difficult problems in terms of structure because depending on their formation history, 
precise composi tion and location, they may b e fluid, solid, or they may even possess a global liquid 



ocean (see iKuchiiMl l2003t iLeger et all , |2004|) 



5 Implications for planetary formation models 

The giant planets in our Solar System have in common possessing a large mass of hydrogen and 
helium, but they are obviously quite different in their aspect and in their internal structures. 
Although studies cannot be conducted with the same level of details, we can safely conclude that 
extrasolar planets show an even greater variety in composition and visible appearance. 

A parallel study of the structures of our giant planets and of giant planets orbiting around 
other stars should provide us with key information regarding planet formation in the next decade 
or so. But, already, some conclusions, some of them robust, others still tentative, can be drawn 
(see also the chapter by Stevenson): 

Giant planets formed in circumstellar disks, before these were completely dissipated: 
This is a relatively obvious consequence of the fact that giant planets are mostly made of hydrogen 
and helium: these elements had to be acquired when they were still present in the disk. Because 
the observed lifetime of gaseous circumstellar disks is of the order of a few million years, this 
implies that these planets formed (i.e. acquired most of their final masses) in a few million years 
also, quite faster than terrestrial planets in the Solar System. 

Giant planets migrated: 

Although not cleanly demonstrated yet, there is evidence that the observed orbital distribution 
of extrasolar planets requires an inward migration of planets, and various mechanisms ha ve been 



proposed for that fsee llda and Linl . l2004al : lAlibert et al.Ll2005l : lMoorhead and A dams. 2005, ...etc.) 



Separately, it was shown that several properties of our Solar System can be explained if Jupiter, 
Saturn, Uranus and Neptune ended up the early formation phase in the presence of a disk with 
quasi-circular orbit, and with Saturn, Uranus and Neptune significantly closer to the Sun tha n 



they are now, and that these three planets subsequently migrated outward ([Tsiganis et al.Ll2005l ). 

Accretion played a key role for giant planet formation: 
Several indications point towards a formation of giant planets that is dominated by accretion 
of heavy elements: First, Jupiter, Saturn, Uranus and Neptune are all significantly enriched in 
heavy elements compared to the Sun. This feature can be reproduced by core-accretion models, 
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for Jupiter and Saturn at least (| Alibert et all 120051 ). Second, the probability to find a giant 



planet arou nd a solar- t ype star (with stellar type F, G or K) is a strongly ris ing function of stellar 
metallicity ( Gonzale3 . 19981 : Santos et al. . 2004; Fischer and Valenti . 20051 ) , a property that i s 

004bl : 



2004bl : lAlibert et all . l2005l ) 



also well- reproduced by standard core accretion models ( Ida and Linl . 
Third, the large masses of heavy elements inferred in some transiting extrasolar planets as well 
as the apparent corr elation be tween mass of hea vy el ements in the plane t and stellar metallicity 
( Guillot et al.l ( 20061 ) ; see also ISato et al. (|2005h and llkoma et alj is a strong indication 

that accretion was possible and that it was furthermore efficient. It is to be noted that none of 
these key pr operties are directly explained b y formation models that assume a direct gravitational 
collapse (see iBosd. 12004 iMaver et all |2004|) . 

Giant planets were enriched in heavy elements by core accretion, planetesimal delivery and/ or 
formation in an enriched protoplanetary disk: 

The giant planets in our Solar System are unambigously enriched in heavy elements compared 
to the Sun, both globally, and when considering their atmosphere. This may also be the case of 
extrasolar planets, although the evidence is still tenuous. The accretion of a central core can explain 
part of the global enrichment, but not that of the atmosphere. The accretion of planetesimals may 
be a possible solution but in the case of Jupiter at least the rapid drop in accretion efhciency as the 
planet reaches appreciable masses (~ 100 or so) implies that such an enrichment would have 
originally concerned only very deep layers, and would requi re a relatively efficie nt upper mixing of 
these elements, and possibly an erosion of the central core (Guillot ei^,|2004'). 

Although not unambiguously explained, the fact that Jupiter is also enriched in noble gases 
compared to the Sun is a key observation to understand some of the processes occuring in the early 
Solar System. Indeed, noble gases are trapped into solids only at very low temperatures, and this 
tells us either that most of the solids that formed Jupiter wer e formed at very low temperature to 



be ab le to trap gases such as Argon, probably as clathrates (jGautier et al. 1,12001 iHersant et al 



2004), or that the planet formed in an enriched disk as it was being evaporated () Guillot and Hueso . 
200d ). 



6 Future prospects 

We have shown that the compositions and structures of giant planets remain very uncertain. This 
is an important problem when attempting to understand and constrain the formation of planets, 
and the origins of the Solar System. However, the parallel study of giant planets in our Solar 
System by space missions such as Galileo and Cassini, and of extrasolar planets by both ground 
based and space programs has led to rapid improvements in the field, with in particular a precise 
determination of the composition of Jupiter's troposphere, and constraints on the compositions of 
a dozen of extrasolar planets. 

Improvements on our knowledge of the giant planets requires a variety of efforts. Fortunately, 
nearly all of these are addressed at least partially by adequate projects in the next few years. The 
efforts that are necessary thus include (but are not limited to): 

• Obtain a better EOS of hydrogen, in particular near the molecular/metallic transition. This 
will be addessed by the construction of powerful lasers such as the NIF in the US and the 
MegaJoule laser in France, and by innovative experiments such as shocks on pre-compressed 
samples. One of the challenges is not only obtaining higher pressures, but mostly lower tem- 
peratures than currently possible with single shocks. The parallel improvement of computing 
facilities should allow more extended numerical experiments. 
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• Calculate hydrogen-helium and hydrogen-water phase diagrams. (Other phase diagrams 
arc desirable too, but of lesser immediate importance). This should be possible with new 

numerical experiments. 

• Have a better yardstick to measure solar and protosolar compositions. This may be addressed 
by the analysis of the Genesis mission samples, or may require another future mission. 

• Improve the values of J4 and Je for Saturn. This will be done as part of the Cassini-Huygens 
mission. This should lead to better constraints, and possibly a determination of whether the 
interior of Saturn rotates as a solid body. 

• Detect new transiting extrasolar planets, and hopefully some that are further from their 
star. The space missions CoRoT (2006) and Kepler (2008) should provide the detection and 
characterization of many tens, possibly hundreds of giant planets. 

• Model the formation and evolution of ice giants such as Uranus, Neptune, and similar planets 
around other stars, in order to analyze detections of these objects and understand planetary 

formation. 

• Improve the measurement of Jupiter's gravity field, and determine the abundance of water in 
the deep atmosphere. This will be done by the Juno mission (launch 2011) with a combination 
of an exquisite determination of the planet's gravity field and of radiometric measurements 

to probe the deep water abundance. 

• It would be highly desirable to send a probe similar to the Galileo probe into Saturn's 
atmosphere. The comparison of the abundance of noble gases would discriminate between 
different models of the enrichment of the giant planets, and the additional measurement of 
key isotopic ratio would provide further tests to understand our origins. 

Clearly, there is a lot of work on the road, but the prospects for a much improved knowledge 
of giant planets and their formation are bright. 
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